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A  STATISTICAL  EVALUATION  OF  SOIL 
AND  CLIMATIC  PARAMETERS  AFFECTING 
THE  CHANGE  IN  PAVEMENT  DEFLECTION 
DURING  THAWING  OF  SUBGRADES 

Edwin  J.  Chamberlain 


INTRODUCTION 

Recent  developments  in  the  design  of  flexible 
pavement  systems  utilize  mechanistic  methods 
rather  than  the  empirical  design  methods  em¬ 
ployed  in  the  past  Rather  than  using  the  results 
of  static  load  tests  such  as  the  plate  bearing  or 
California  bearing  ratio  (CBR)  tests,  the  new  de¬ 
sign  methods  require  the  input  of  resilient  modu¬ 
li  for  each  of  the  layers  in  a  pavement  system 
While  preparing  a  recommended  procedure  for 
determining  the  resilient  moduli  of  base,  sub¬ 
base,  and  subgrade  materials  during  thawing, 
the  author  reviewed  literature  for  both  labora¬ 
tory  and  field  studies. 

One  field  study  that  was  of  special  promin¬ 
ence  was  reported  by  Scrivner  et  al  (1969)  Be¬ 
cause  this  study  provided  a  considerable  amount 
of  data  for  a  variety  of  soil  and  climatic  condi¬ 
tions,  several  of  the  potentially  significant  par¬ 
ameters  affecting  the  change  in  pavement  per¬ 
formance  could  be  evaluated  Two  factors  were 
of  interest:  1)  the  increased  resilient  pavement 
deflection  during  thaw,  and  2)  the  time  required 
for  recovery  of  the  resilient  pavement  deflection 
to  prefrozen  levels  In  order  to  evaluate  the  in¬ 
fluence  of  soil  and  climatic  parameters  on  these 
factors,  the  author  employed  statistical  tech¬ 
niques  The  results  were  somewhat  surprising 
<md  are  presented  here  to  stimulate  discussion 

Scrivner  et  al  (1969)  applied  a  cyclic  load  to 
existing  flexible  pavement  surfaces  using  the  Dy- 
naflect,  a  pavement-mounted,  cyclic  pavement¬ 
loading  device  that  applies  a  load  of  1,000  lb  at 
a  frequency  of  8  Hertz,  and  recorded  the  recov¬ 
erable  or  resilient  surface  deflection  Six  test 
sites  at  each  of  four  areas  between  Springfield, 
Illinois,  and  Duluth,  Minnesota,  were  observed 
over  a  1-year  period  The  mean  freezing  index 


was  100°F-days  in  the  most  southerly  of  the 
regions  and  2100°F-days  in  the  most  northerly. 
Where  the  freezing  index  was  small  (100°F-day$) 
the  frost  penetration  into  the  subgrade  was  shal¬ 
low  and  intermittent,  while  where  the  freezing 
index  was  large  (2100°F-days),  the  subgrade 
freezing  was  deep  and  uninterrupted.  Figure  1 
gives  typical  frost  penetration  for  each  climatic 
region. 

This  study  analyzes  the  influence  of  the  freez¬ 
ing  index,  the  number  of  freeze-thaw  cycles,  the 
depth  of  frost  penetration  into  the  subgrade,  and 
the  subgrade  soil  type  on  the  seasonal  change  of 
resilient  deflection  due  to  thawing  and  on  the 
time  after  thaw  to  recovery  of  the  normal  pave¬ 
ment  deflection 


DESCRIPTION  OF  TEST  SITE 

As  shown  in  Table  1,  a  wide  range  of  subgrade 
materials  was  encountered,  from  clay  tills  to 
sand  fills,  with  silty  clays  being  the  most  com¬ 
mon  Base  and  subbase  materials  included 
native  granular  material,  crushed  stone,  gravel, 
crushed  rocks,  sand-gravel,  and  NaCI-stabilized 
and  bituminous-treated  gravels.  Pavements  were 
predominantly  asphalt  concrete  A  few  surface 
treatment  pavements  were  also  encountered 

In  each  of  the  four  study  areas  two  of  the  six 
test  sections  were  located  on  highways  designed 
to  carry  relatively  heavy  traffic,  and  four  were 
located  on  roads  designed  for  light  traffic  The 
test  sections  were  1000  feet  in  length.  Pavement 
deflections  at  ten  test  points  in  the  outer  wheel 
path  were  observed  over  a  period  from  Decem¬ 
ber  1966  to  October  1967  Approximately  20 
visits  were  made  to  each  test  section  during  this 
period 
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Pavement  deflection  and  frost  penetration  vs  time.  Typical  plots  tor  freezing  indexes  of  a)  100 ,  b)600f  c)  1300 ,  and 
} F-days  (from  Scrivner  et  al.  1 969). 


Table  1.  Properties  of  test  sections  (from  Scrivner  et  al.  1969). 


ARI  A 

SI  (  . 

NO. 

m  SKiN 

V  1  ASS.  ' 

IliU  KM  SS  (INCHIS) 

MAH  RIA1 

TVl’l 

SI  Rl  AC'} 

HASI 

si  bbasi: 

SI Rl  AC  1 

BAS! 

SL  BBASI 

SI  BORADI 

l 

19 

1 

1  *' 

7 

_ 

ST. 

Ciran.  mall. 

_ 

Silty  clay 

1 

20 

1 

4.5 

8 

b 

A.C. 

Cr.  stone 

Cr.  stone 

Silty  clay 

1 

21 

1 

2.5  '* 

6.5 

— 

ST. 

Cr.  stone 

— 

Silty  clay 

t 

i2 

1 

4.3  b 

b 

— 

ST. 

Ci  ravel 

— 

Silty  clay 

1 

23 

1 

i  «• 

7 

— 

ST. 

Cr.  stone 

— 

Silty  clay 

1 

24 

2 

4.5 

9 

b 

A.C 

Cr.  stone 

Ciravel 

Silty  clay 

i 

13 

1 

0.5  '■ 

8 

— 

ST. 

Ci  ravel 

— 

Silty  clay 

2 

14 

1 

0.5  " 

4 

4 

SI. 

NaCI  stab.  gr. 

Ciravel 

Silty  clay 

i 

15 

1 

2 

8 

— 

A.C. 

Gravel 

— 

Clay  (ill  (some  gravel ) 

2 

lb 

1 

0.5  1 

4 

4 

SI. 

Ci  ravel 

Ciravel 

Silty  clay 

i 

17 

1 

3 

10 

— 

A.C. 

Cr.  stone 

— 

Sand  fill 

2 

18 

1 

2 '' 

7" 

— 

A.C. 

Ci  ravel 

— 

Silty  clay 

■ 

1 

1 

3 

3 

9 

A.C. 

Cr.  rock 

S.tnd-gravel 

Silty  day  loam 

3 

2 

I 

4 

lb.  5 

— 

A.C. 

Gravel 

— 

Silty  clay  loam 

3 

3 

1 

2 

3 

9 

A.C. 

Cr.  rock 

Sand-gravel 

Organic  loam 

3 

4 

1 

1.5 

9 

— 

A.C. 

Gravel 

— 

Organic  loam 

3 

5 

2 

3 

b 

12 

A.C 

Cr.  rock 

Sand-gravel 

Silly  loam 

3 

b 

i 

b 

4.5 

12 

A.C. 

Gravel 

Sand-gravel 

Clay  till  (some  gravel ) 

4 

7 

1 

b 

13 

_ 

A.C 

Gravel 

_ 

A-7-5  clay 

4 

8 

i 

6 

lb.5 

— 

A.C. 

Gravel 

— 

Gravel  fill  on  sand>  loam 

4 

9 

1 

4 

3 

9 

A.C. 

Ciravel 

Ciravel 

Gray  and  red  day 

4 

10 

I 

4 

9 

— 

A.C. 

Sand-gravel 

— 

Clay  loam 

4 

1  1 

1 

3 

3 

9 

A.C. 

Ciravel 

Sand-grave) 

A-b  clay  to  A-4  sandy  loam 

4 

12 

i 

8 

4 

12 

AC 

Hit.  treat,  gr 

Ci  rave1 

Sand  fill 

*  l 

-Design 

for  light  traffic 

2  Design  to 

i.irrs  rcl.it i 

vvly  hf.ivy 

traffic. 

h  Thickness  determined  by  drilling 


Thermocouples  were  installed  at  two  loca¬ 
tions  within  each  test  section  The  depth  of  frost 
penetration  was  determined  from  the  position  of 
the  32°  isotherms  and  averaged  for  each  of  the 
two  locations  in  each  test  section 


TEST  RESULTS 

Figure  1  shows  typical  results  for  each  of  the 
four  study  areas.  Figure  la  shows  that  in  the 
most  southerly  locations  (freezing  index  = 
100°F-days),  the  frost  penetration  into  the  sub¬ 
grade  was  shallow  and  intermittent.  Several 
freeze-thaw  cycles  occurred  over  the  winter 
period.  Farther  north,  where  the  freezing  index 
was  600°F-days,  the  frost  penetration  was  still  in¬ 
termittent  but  the  frost  depth  was  significantly 
greater.  Still  farther  north  (1300°F-days  freezing 
index),  the  frost  penetration  was  more  or  less 
continuous  over  the  winter  and  approached  40 
inches  At  the  most  northerly  location  (2100°F- 
days  freezing  index),  the  frost  penetration  into 
the  subgrade  was  uninterrupted  and  exceeded 
60  inches 

In  each  of  these  representative  plots  the 
spring  thaw  is  accompanied  by  an  upward  spike 
in  the  Dynaflect  deflection  curve  Generally,  the 


ratio  of  the  magnitude  of  this  spike  to  the  nor¬ 
mal  fall  deflection  value  appeared  to  be  unre¬ 
lated  to  freezing  index  in  the  examples  shown  in 
Figure  1  Only  the  most  northerly  site  is  an  ex¬ 
ception  to  this  observation. 

Figure  1  shows  that  the  time  to  recovery  was 
relatively  short  The  area  under  the  spike  gener¬ 
ally  covered  a  period  of  approximately  one 
month.  The  spike  is  generally  followed  by  a 
slight  increase  in  deflection  in  the  next  month 
and  finally  a  gradual  decrease  over  the  re¬ 
mainder  of  the  summer 

Table  2  presents  the  test  results  for  all  test 
sites  The  first  column  gives  the  section  identifi¬ 
cation  number  according  to  Scrivner  et  al.  (1969). 
The  next  column  gives  the  freezing  index  (/),  also 
as  reported  by  Scrivner  et  al  The  third  and 
fourth  columns  give  the  number  of  freeze-thaw 
cycles  (N)  and  the  depth  of  frost  penetration  (D) 
into  the  subgrade  as  determined  by  the  author 
from  the  reported  data  plots.  The  next  column 
gives  my  estimate  of  the  frost  classification  (F)  of 
the  subgrade  material  according  to  the  Corps  of 
Engineers  system  (U  S  Army  1965)  where  1  is  the 
lowest  frost  susceptibility  class  and  4  is  the  high¬ 
est  frost  susceptibility  class. 

The  next  two  columns  show  the  pavement  de¬ 
flection  data  reported  by  Scrivner  et  al  (1969) 
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The  normal  fall  deflection  data  (W1n)  were  taken 
directly  from  tabulated  data  in  their  report  while 
the  maximum  spring  deflection  (W1max)  was  in¬ 
terpreted  from  data  plotted  in  their  report.  They 
define  normal  deflection  as  "the  average  deflec¬ 
tion  observed  during  the  months  of  August,  Sep¬ 
tember,  and  October  prior  to  the  first  freeze  " 
Both  Wln  and  VVlmdx  are  average  values  deter¬ 
mined  from  10  locations  within  each  test  sec¬ 
tion  The  next  column  shows  the  percent  in¬ 
crease  in  deflections  (A)  during  the  spring  period 
The  last  column  shows  the  time  (t2 0)  from  the 
onset  of  thawing  to  recovery  of  all  but  20%  of 
Wln  This  recovery  time  was  arbitrarily  selected 
to  eliminate  the  long,  slow  recovery  period  lead¬ 
ing  into  the  fail  season,  which  is  of  minor  conse¬ 
quence  in  the  performance  of  the  pavement 
after  th^w 


METHOD  OF  ANALYSIS 

My  analysis  assumes  that  the  increase  in  pave¬ 
ment  deflection  is  wholly  attributable  to  the 
weakening  of  the  subgrade  as  a  result  of  freezing 
and  thawing  Obviously,  changes  in  the  stiffness 
of  the  pavement  and  the  base  and  subbase  mate¬ 
rial  also  probably  occurred  However,  neither  of 
these  factors  could  be  quantified  for  this 
analysis,  so  they  are  not  directly  considered 

Data  were  also  unavailable  on  several  other 
potentially  significant  factors,  including  mois¬ 
ture  content,  dry  density,  proximity  of  the  water 
table,  and  precipitation  Fven  so,  it  was  felt  that 
any  sensitivity  to  the  severity  of  the  climate  and 
the  other  factors  available  could  be  identified 
To  do  this,  empirical  relationships  between  the 
dependent  variables  A  and  f,0  and  the  indepen¬ 
dent  variables  Wln.  I,  N.  D,  and  F  were  assumed 
(A  was  also  considered  as  an  independent 
variable  tn  evaluating  f20)  The  validity  of  these 
relationships  was  tested  by  multiple  linear  re¬ 
gression  analysis  and  analysis  ot  variance 

A  visual  inspection  of  the  data  revealed  that 
the  independent  variables  might  be  linearly 
related  to  the  dependent  variables  Therefore, 
the  following  expressions  were  tried 

A  =  a0  +  a,Wln  +  a2/  +  a  }N  +  aAD  +  a^F  (1) 

and 


where  the  subscripted  terms  a  and  b  are  the  re¬ 
gression  coefficients  Second  degree  governing 
equations  were  also  analyzed  to  allow  for  non¬ 
linear  relationships,  but  the  results  were  not  sig¬ 
nificantly  different  Only  the  results  of  the  first 
degree  equations  are  reported  here 

The  statistical  analysis  employed  has  been  de¬ 
scribed  by  Robbins  and  Chamberlain  (1973)  It 
employs  a  multiple  linear  regression  analysis  of 
the  functional  relationships  between  the  inde¬ 
pendent  and  dependent  variables  of  interest 
while  allowing  the  successive  selection  of  the 
most  influential  independent  variables  with  re 
spect  to  the  variance  reduction  of  the  dependent 
variable  The  acceptance  or  rejection  of  a  vari¬ 
able  was  governed  by  an  F  test  in  which  the  max¬ 
imum  probability  level  (a)  for  rejecting  a  vari¬ 
able  could  be  selected  as  desired 

A  5%  a  level  was  used  The  a  value  governing 
the  entering  or  removal  of  an  independent  vari 
able  allowed  the  testing  of  the  contribution  of 
that  variable  to  the  reduction  in  variance  of  the 
dependent  variable  in  the  presence  of  any  other 
variables  previously  accepted  Once  a  new  vari¬ 
able  was  accepted,  all  other  accepted  variables 
were  subjected  to  the  F  test  and  any  that  were 
rejected  were  removed  and  placed  in  the  list  to 
be  reconsidered  The  process  of  judging  accep¬ 
tance  and  rejection  continued  until  no  more 
variables  passed  the  F  test 

The  statistical  analysis  also  provided  a  meas 
ure  of  the  strength  of  the  linear  relationship  be¬ 
tween  variables  through  the  correlation  coeffi¬ 
cient  r  If  the  data  points  lay  on  a  straight  line, 
the  correlation  coefficient  would  be  ±  1.  If  the 
data  points  were  randomly  scattered,  the  corre¬ 
lation  coefficient  would  be  0  Another  statistical 
element,  the  standard  error,  provided  a  measure 
of  the  scatter  of  the  data  about  the  regression 
line  For  normally  distributed  data,  approxi¬ 
mately  two-thirds  of  the  data  points  fall  within  ± 
one  standard  error,  while  95%  fall  within  ±  two 
standard  errors 


RESULTS  OF  ANALYSIS 

Change  in  resilient  deflection 
due  to  thawing,  A 

The  results  of  the  regression  analysis  (Table  3) 
for  A  revealed  that  only  the  frost  classif ication  F 
was  significant  at  a  probability  of  occurrence 
level  (1-a)  of  95%  or  better  The  correlation  co¬ 
efficient  was  0  50  and  the  standard  error  39  55 
Figure  2a  illustrates  the  data  of  Scrivner  et  al 


(2) 


Table  3.  Results  of  regression  analysis. 


Independent 

Probability  oi 

Independent 

variables 
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Equation  coefficients 

Dependent 

variables 

accepted 

a  variable 

Intercept  Slope 

Correlation  Standard 

variable 

entered 

o  tt  5  % 

o(%) 

a(> 

coefficient  error 

A% 

Wln.  I.  N .  D,  F 

F 

1 

-42  69% 

32  905%/unit 

0  50 

39  55 

A 

w1n 

not  passing 

31 

39  00% 
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44  52 

A 
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not  passing 

52 
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013 

45  20 

A 

N 

not  passing 

>75 

61  68% 

1  063%/umt 
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A 
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F 
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Wln.  /,  N,  D,  F,  A 

D 

1 

52  59  days 

-0  144  day/inch 

0  51 

5  10 

*20 

wln 

yvln 

5 

42  29  days 

3  450  days/milli-inch 

0  41 

5  42 

*20 

/ 

not  passing 

6 

51  48  days 

-0  00299  day/°F-day 

040 

5  45 

*20 

N 

not  passing 

7 

45  95  days 

1  091  days/unit 

0  38 

5  50 

*20 

D 

D 

1 

52  59  days 
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0  51 

5  10 

*20 

F 

not  passing 

37 

43  00  days 

1  667  days/unit 

019 

5  83 

*20 

A 

not  passing 

30 

46  56  days 

0  0288  day/% 

0  22 

5  79 

Figure  2  Data  points  and  regression  line  for  A 
and  t20  versus  F 


(1969)  and  the  best  fit  line  in  the  A  versus  /  plane 
There  is  a  sharp  increase  in  A  with  increasing 
frost  susceptibility  and  considerable  scatter  in 
the  data,  reflected  by  the  relatively  low  correla¬ 
tion  coefficient  and  the  high  standard  error. 

No  other  variable  passed  the  desired  F  test  of 
a  =  5%  In  fact,  none  of  the  remaining  indepen¬ 
dent  variables  passed  an  F  test  for  a  =  25%  In 
the  presence  of  all  other  parameters,  the  most 
influential  parameter  on  A  is  the  frost  classifica¬ 
tion;  the  inclusion  of  any  of  the  other  variables 
would  not  contribute  significantly  to  the  de¬ 
creasing  of  the  variance  in  A 

It  was  also  desired  to  test  the  individual  con¬ 
tributions  of  each  independent  variable  on  A 
The  functional  relationships  were  assumed  to  be 
of  the  form 

A  =  c0  +  ClX  (3) 

where  X  is  the  independent  variable  of  interest 
and  cQ  and  c1  are  the  regression  coefficients  The 
results  of  this  analysis  are  given  in  Table  3  The 
frost  susceptibility  classification  of  the  subgrade 
remained  the  most  significant  parameter  The 
depth  of  frost  penetration  was  the  least  signifi¬ 
cant.  The  probability  of  any  variable  other  than 
the  frost  classification  significantly  influencing 
the  reduction  in  variance  in  A  was  no  better  than 
30%,  so  all  other  variables  were  rejected  by  the 
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Figure  3.  Data  points  and  regression  fine  for  A  and 
t2n  versus  W/n. 
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Figure  4  Data  points  and  regression  line  for  A  and 
t  ,0  versus  I 


F  test 

The  linear  fits  established  by  the  regression 
analysis  are  plotted  along  with  the  data  of  Scriv- 
ner  et  al  (1969)  in  Figures  2a-6a  The  following 
trends  can  be  observed: 

1  A  increases  significantly  with  increasing  F 

2  A  increases  with  increasing  Wln 

i  A  increases  slightly  with  increasing  / 

4  A  increases  slightly  with  increasing  N 

5  A  increases  slightly  with  increasing  D 
However,  since  the  probabilities  of  falsely  ac¬ 
cepting  the  variables  W1n,  I,  N,  and  D  are  so 
great,  the  use  of  their  linear  fits  is  questionable 

Recovery  time  after  onset  of  thawing,  f20 

The  results  of  the  statistical  analysis  for  f2Q 
(Table  3)  show  that  the  depth  of  frost  penetra¬ 
tion  D  is  the  most  significant  parameter  in  the 
presence  of  the  others  (a  =  1%)  and  that  no 
other  variable  passes  the  F  test  for  a  =  5% 
When  considered  separately,  the  variables  W1n, 
/,  and  N  pass  F  tests  a  =  5,  6,  and  7%,  respec¬ 
tively  The  a  levels  for  the  remaining  two  par 
ameters,  A  and  the  frost  classification,  were  very 
high,  so  they  are  not  considered  statistically  sig¬ 
nificant 


As  illustrated  in  Figure  6b,  the  recovery  time 
f20  decreases  with  increasing  depth  of  frost  pene¬ 
tration  D  The  fit  of  D  to  f20  appears  to  be  a  little 
better  than  that  of  F  to  A  The  correlation  coeffi¬ 
cient  is  0  51  and  the  standard  error  is  5  10  The 
results  of  the  individual  linear  fits  between  f20 
and  the  independent  variables  are  given  in  Table 
3  and  illustrated  in  Figures  2b-6b  and  7  along 
with  the  data  of  Scrivner  et  al.  (1969)  The  follow¬ 
ing  trends  were  observed: 

1  t2()  increases  slightly  with  increasing  F 

2  t20  increases  slightly  with  increasing  W1n 

3  f20  decreases  slightly  with  increasing  / 

4  f20  Increases  slightly  with  increasing  N 

5  f2()  decreases  with  increasing  D 

6  f20  increases  slightly  with  increasing  A 

DISCUSSION  OF  RESULTS 

Change  in  resilient  deflection 
due  to  thawing,  A 

Of  special  interest  in  this  analysis  is  the  obser¬ 
vation  that  the  percent  increase  in  pavement  de¬ 
flection  upon  thawing  is  not  significantly  related 
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Figure  5.  Data  points  and  regression  line  for  A  and 
versus  N. 
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Figure  7.  Data  points  and  regression  line  for  ijo 
versus  A 
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Figure  6.  Data  points  and  regression  line  for  A  and 
t jo  versus  D 

to  the  freezing  index.  Since  a  similar  observation 
was  made  by  the  Canadian  Good  Roads  Associa¬ 
tion  (1965)  from  the  results  of  Benkelman  beam 
rebound  tests  in  Canada,  it  can  be  assumed  that 
this  result  of  the  analysis  is  reasonable  The  sig¬ 
nificance  of  this  observation  is  that  if  frost  pene¬ 
trates  into  the  subgrade  of  F 2,  F3,  or  FA  frost 
classification,  7  able  2,  the  potential  loss  in  pave¬ 
ment  carrying  capacity  is  the  same  regardless  of 
the  seventy  of  the  freezing  index 

The  ranking  by  o-level  for  Wln,  /,  N  and  D  re¬ 
vealed  that  o  =  1,  31,  52,  and  >75%,  respective¬ 
ly  Because  the  probabilities  of  rejecting  all  vari¬ 
ables  but  F ,  even  when  variables  were  con¬ 
sidered  separately,  were  30%  or  better,  only  F 
can  be  considered  as  a  predictor  of  A  The  reli¬ 
ability  of  using  F  as  a  predictor  for  A  is  poor  as 
revealed  by  the  large  standard  error  The  scatter 
may  be  a  result  of  the  interpretation  required  to 
extract  the  data  from  the  report  of  Scrivner  et  al 
(1969)  or,  more  likely,  the  result  of  the  unavaila¬ 
bility  of  all  significant  variables  for  the  analysis. 
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That  the  most  significant  variable  influencing 
the  reduction  in  variance  in  A  was  the  frost  clas¬ 
sification  confirms  the  need  to  carefully  select 
the  soil  materials  for  pavement  structures  and  to 
provide  a  means  of  accurately  identifying  the 
frost  weakening  characteristics  of  soils.  Because 
the  weakening  of  asphaltic  pavement  structures 
during  thawing  is  a  system  response  involving 
both  surfacing,  and  compacted  and  natural  soil 
materials,  any  reliable  frost-weakening  criteria 
must  be  developed  in  terms  of  the  in-place  con¬ 
ditions  of  each  layer 

The  most  direct  approach  to  this  analysis 
would  require  in-place  measurement  of  pave¬ 
ment  deflections  under  repeated  loading  along 
with  measurement  of  the  above-mentioned  prop¬ 
erties,  .and  laboratory  repeated-load  tests  on 
each  of  the  components  of  the  pavement  sys¬ 
tem  A  layered  system  mathematical  analysis 
would  be  required  to  relate  the  laboratory  obser¬ 
vations  to  the  field  performance  These  methods 
are  presently  being  used  at  CRREL  (Johnson  et 
al.  1978)  If  a  sufficient  number  of  laboratory 
and  field  observations  were  correlated,  the  labo¬ 
ratory  and  field  test  could  probably  be  replaced 
or  supplemented  by  common  soil  index  tests 
such  as  grain  size,  Atterberg  limits,  dry  density, 
moisture  content,  permeability,  and  soil  suction 

Recovery  time  after  thawing,  t20 

An  unanticipated  result  of  this  analysis  is  that 
the  recovery  time  decreases  with  depth  of  freez¬ 
ing  and  the  freezing  index;  that  is,  the  greater  the 
freezing  index,  the  greater  will  be  the  depth  of 
freezing  and  the  shorter  the  recovery  time.  This 
observation  originally  conflicted  with  the 
author's  intuitive  judgment  that  the  recovery 
time  would  be  longer  with  increasing  frost  depth 
because  the  greater  frost  depth  would  require  a 
longer  thaw  period  and  would  cause  a  change  in 
deformational  properties  over  a  greater  depth  A 
closer  look  at  the  frost  penetration- time  curves 
(Fig  1 )  and  recovery  time  data  in  T able  2  reveals 
that  the  thaw  period  (a  few  days  to  two  weeks)  is 
relatively  short  in  comparison  with  the  recovery 
period  (35  to  60  days)  and,  thus,  is  not  the  domi¬ 
nant  factor  in  determining  the  length  of  the  re¬ 
covery  period 

A  possible  explanation  for  the  shorter  re¬ 
covery  period  occurring  for  the  case  of  deepest 
frost  penetration  is  that  the  freezing  and  thawing 
process  changes  the  particle  arrangement  in  the 
subgrade  soil  and  gives  it  a  higher  permeability, 
thus  improving  the  drainage  Recent  tests  by  the 
author  (Chamberlain  and  Blouin  1978,  Chamber- 


Figure  8.  Permeability  vs  void  ratio  for  Toledo 
Penn  7  dredged  material.  Numbers  adjacent  to 
points  indicate  the  number  of  freeze-thaw 
cycles.  (From  Chamberlain  and  Cow  1978) 


lain  and  Gow  1978)  revealed  that  the  permeabili¬ 
ty  of  a  soil  can  be  increased  by  as  much  as  two 
orders  of  magnitude  by  frost  action  Figure  8,  for 
example,  shows  the  large  increase  in  permeabil¬ 
ity  for  a  dredged  material  It  appears  that  the 
shorter  recovery  time  for  the  greater  depth  of 
freezing  may  be  a  consequence  of  the  improved 
drainage  conditions  imparted  by  frost  action 
The  frost  penetration  depth  D  was  by  far  the 
most  significant  variable  affecting  t20  in  the  pres¬ 
ence  of  all  others  The  probability  of  falsely  ac¬ 
cepting  this  variable  is  1%.  In  the  presence  of  D 
no  other  variable  would  pass  a  significance  level 
of  even  a  =  20%,  so  all  other  variables  were  re¬ 
jected  However,  when  considered  singly,  the 
Wln,  /,  and  N  variables  passed  the  5%,  6%,  and 
7%  significance  levels  respectively  This  means 
that  W1n,  I,  and  N  correlate  strongly  with  D  and 
are  good  predictors  of  f20,  that  given  the  choice 
of  these  four  variables  D  would  be  the  best 
predictor  of  f20,  and  that  including  W1n,  /,  and/or 
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N  would  not  contribute  significantly  to  the 
prediction  of  f20  in  the  presence  of  D;  i.e.,  the  ac¬ 
curacy  of  the  equation 

*20  =  a0  +  a1D 

in  predicting  f20  is  not  improved  significantly  by 
including  the  variables  W1n,  /,  or  N. 

Even  when  considered  singly  in  the  statistical 
analysis,  the  variables  F  and  A  fail  to  pass  a  sig¬ 
nificant  level  of  25%,  and  thus  cannot  be  con¬ 
sidered  as  reliable  predictors  or  indicators  of  t2o 

The  reliability  of  D  as  a  predictor  of  f20  is  good 
as  revealed  by  the  relatively  small  standard  error 
on  the  regression  line  between  D  and  t20.  How¬ 
ever,  because  of  the  limited  number  of  sites  eval¬ 
uated,  the  results  should  not  be  interpreted  as 
establishing  a  universal  relationship  between  D 
and  f20.  The  analysis  should  be  tested  on  a  great 
many  more  sites 


CONCLUSIONS 

The  most  interesting  observations  of  this  sta¬ 
tistical  analysis  are  that  the  freezing  index  is  not 
a  significant  parameter  in  determining  the  percent 
increase  in  pavement  deflection  during  thawing 
and  that  the  recovery  time  decreases  with  the 
depth  of  freezing. 

The  necessity  for  careful  selection  of  soil 
materials  for  pavement  systems  is  reflected  in 
the  observation  that  the  most  significant  varia¬ 
ble  affecting  the  reduction  in  variance  in  A  was 
the  frost  susceptibility  classification 

Since  the  standard  error  was  large  for  the 


linear  fit  between  A  and  F,  the  established  rela¬ 
tionship  may  not  give  an  accurate  estimate  of  A 
from  knowledge  of  F.  Much  more  study  is  re¬ 
quired  to  analyze  this  relationship  and  other 
variables  such  as  moisture  content,  dry  density, 
and  depth  to  water  table  must  be  included 
Although  the  recovery  time  was  found  to  be 
strongly  related  to  the  depth  of  freezing,  the 
number  of  sites  evaluated  was  not  large  A  great 
many  more  sites  must  be  evaluated  to  provide  a 
full  test  of  the  results 
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